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f— I . Abstract. In this paper we present some test results of our newly developed 

' Multi-Phase Chemo-Dynamical Smoothed Particle Hydrodynamics (MP-CD-SPH) 

code for galaxy evolution. At first, we present a test of the "pure" hydro SPH 

■ part of the code. Then we describe and test the multi-phase description of the 
' gaseous components of the interstellar matter. In this second part we also compare 

our condensation and evaporation description with the results of a previous 2d 
multi-phase hydrodynamic mesh code. 
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Since several years Smoothed Particle Hydrodynamics (SPH) (Mon- 

■ aghan, 1992) has been applied successfully to study the formation and 
evolution of galaxies. Its Lagrangian nature as well as its easy imple- 
mentation together with standard N-body codes allows for a simul- 
taneous description of complex dark matter-gas-stellar compositions 
(Navarro & White, 1993; Mihos & Hernquist, 1996). Nevertheless, un- 
til now the codes lack processes that are based on the coexistence of 
different phases of the interstellar medium (ISM), mainly dissipative, 
dynamical and stellar feedback, element distributions, etc. We have 
therefore developed a 3d chemo-dynamical code which is based on 
our single phase galactic evolutionary program (Berczik, 1999; Berczik, 
2000). This code includes many complex effects such as a multi-phase 
ISM, cloud-cloud collisions, the drag force between different ISM com- 
ponents, condensation and evaporation of clouds (CE), star formation 
(SF) and stellar feedback (FB). The more detailed description of the 
new code will be presented in a comprehensive paper (Berczik et al., 
in preparation). 
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2. The hydro code test 



The self-gravitating collapse of an initially isothermal, cool gas sphere 
is a common test problem for SPH codes (Evrard, 1988; Hernquist & 
Katz, 1989; Steinmetz & Miiller, 1993; Carraro et al., 1998; Thacker et 
al., 2000; Springel et al., 2001). Following these authors, we consider a 
gas sphere of total mass M, radius R, and initial density profile p(r) = 
2 f and with an internal energy per unit mass of u = 0.05 ^-jf-- 
At the beginning, the gas particles are at rest. We use a system of units 
with G=M=R=1. In Fig. 1. the time evolution of the different energies 
(left) and of the relative error of the total energy (right) are displayed. 
In the right figure we also compare results applying the publicly accessi- 
ble GADGET code (Springel et al., 2001) with "standard" parameters 
for the 32,000 particles. During the central bounce around t « 1.1 most 
of the kinetic energy is converted into heat and a strong shock wave 
travels outwards. The results of Fig. 1. agree very well with those of 
Steinmetz & Miiller (1993) and Springel et al. (2001). The maximum 
relative total energy error is around 0.05 % even for low (8,000) particle 
numbers. 




Figure 1. Time evolution of the thermal, kinetic, potential and total energy for the 
collapse of an initially isothermal gas sphere (left) and the relative total energy error 
for different models (right). The different lines correspond to the different numbers 
of gas particles. In the right figure we also present the energy error of the GADGET 
code with "standard" parameters for the 32,000 particles. 



3. The multi-phase processes 

In our multi-phase gas code we use a two-component gas description of 
the ISM (Theis et al., 1992; Samland et al., 1997), a cold (10 2 -10 4 K) 
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cloudy component treated by sticky particles (Theis & Hensler, 1993) 
and a diffuse hot gas (10 4 -10 7 K) represented by SPH particles. The cold 
clumps follow an empirical mass-radius relation (Larson, 1981; Solomon 
et al., 1987; Maloney, 1 990; Inoue k K amaya, 2000; Theis & Hensler, 
1993) with h cl ~ 50 v /m cl /(10 ti M ) pc and experience dissipation 
due to cloud-cloud collisions and drag forces within the hot gas. The 
mass exchange between "cold" and "hot" gas phases happens basically 
through condensation and/or evaporation (CE) described by Cowie et 
al. (1981) and Koppen et al. (1998). The basic parameter switching 
between both processes, ao, is set to 0.03 according to Cowie et al. 
(1981). 




t (Gyr) t (Gyr) 

Figure 2. The time evolution of the mass exchange rate due to condensation vs. 
evaporation (CE) effects (left). The letters indicate the phase (i.e. H - Hot; C - 
Cold) and the numbers indicate the different models. In the right figure we compare 
our "best" model with 2d multi phase hydro mesh code results, for similar initial 
conditions. 

As a CE test for our code, we calculate the evolution of a "hypotheti- 
cal" gas system inside a dark matter (DM) halo and compare the results 
with a 2d multi-phase hydro mesh-code (Rieschick & Hensler, 2000) for 
the same system. The initial total gas mass amounts to 2 x 10 9 M 
(99% "COLD" + 1% "HOT") and is distributed according to a Plummer- 
Kuzmin disk with parameters a=0.1 kpc and b=2 kpc (Miyamoto & 
Nagai, 1975). The initial temperature for the cold gas is set to 2xl0 3 
K, for the hot gas to 10 6 K. The DM halo is static with parameters 
ro=2 kpc and po = 2 M /pc 3 (Burkert, 1995) so that its mass inside 
the initial gas distribution (20 kpc) amounts to ~ 3.2 x 10 11 M . 

In some of the models we also use an initial cloud mass function 
(ICMF) (Miyazaki & Tsuboi, 2000) with a power-law index -1.5 and 
lower and upper mass limits at 10 4 M and 10 6 M , respectively. In 
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Fig. 2. we present the time evolution of the mass exchange rate via 
CE for different models: VER-01 with 5,000 hot-gas SPH particles of 
identical mass and 5,000 cold gas particles, VER-02 with 10,000 particles 
each, VER-03 with the same particle numbers as in VER-01 but an 
ICMF, and VER-04 equal to VER-02 but also with an ICMF. The figures 
show that the basic behaviour of mass exchange between the hot and 
cold phase is well described with our model and almost independent of 
the particle number and ICMF, but in agreement on average with the 
rates of the 2d mesh code. 
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